Low-temperature phases in PbZro.52Tio.48O3 : A neutron powder diffraction study 
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A neutron powder diffraction study has been carried out on PbZro.52Tio.48O3 in order to resolve an 
ongoing controversy about the nature of the low-temperature structure of this strongly-piezoelectric 
and technologically-important material. The results of a detailed and systematic Rietveld analysis 
at 20 K are consistent with the coexistence of two monoclinic phases having space groups Cm and 
Ic respectively, in the approximate ratio 4:1, and thus support the findings of a recent electron 
diffraction study by Noheda et al. [Phys. Rev. B 66, 060103 (2002)]. The results are compared to 
those of two recent confficting neutron powder diffraction studies of materials of the same nominal 
composition by Hatch et al. [Phys. Rev. B 65, 212101 (2002)] and Frantti et al. [Phys. Rev. B 
66, 064108 (2002)]. 

PACS numbers: 61.12.Ld, 61.50.Ks, 61.66.Fn, 77.84.Bw 
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I. INTRODUCTION 

The strongly piezoelectric system PbZri-j^Tia^Oa 
(PZT) has long been known to have a perovskite-type 
structure with regions of rhombohedral and tetragonal 
symmetry below the ferroelectric Curie point separated 
by an almost vertical line at x « 0.5 in the temperature- 
composition phase diagram, which is known as the mor- 
photropic phase boundary (MPB),^ Following the recent 
discovery of a narrow region with monoclinic Cm symme- 
try in the vicinity of the MPB,^ numerous experimental 
and theoretical studies of PZT and related systems have 
been undertaken in order to clarify the relationships be- 
tween the structural features and the piezoelectric prop- 
erties. As previously discussed;^ the ferroelectric polar- 
ization in the new phase is no longer constrained by sym- 
metry to lie along a symmetry axis, but instead is free to 
rotate within the symmetry plane. Furthermore, because 
of the near-degeneracy of the free energies of the various 
phases, rotation of the polarization axis away from the 
polar axes of the rhombohedral and tetragonal phases can 
be accomplished with an applied electric field, resulting 
in an induced monoclinic phase and a large electrome- 
chanical responsei^ 

The phase diagram of the PZT system around the 
MPB as reported in a recent paper by Noheda et al.^ 
is shown in Fig. 1. Above the Curie temperature, the 
structure is cubic over the entire range of composition, 
with space group Pm3m and lattice parameter ao ~ 4 
A. The rhombohedral region is characterized by high- 
and low-temperature phases (Rht and Rlt) in which 
there are polar shifts of the atoms along the pseudocu- 
bic [111] axis.-'^ R-ht has space group symmetry i?3m, 
with lattice parameters gr ~ oq, and a slightly less 
than 90° (hexagonal values an ~ clq^/2,ch ~ ao%/3). 
In Rlt there are additional displacements of the oxy- 
gen atoms superimposed on the ferroelectric shifts due 
to antiphase tilting of the oxygen octahedra about the 
[111] axis, corresponding to an R-point instability. As 
a consequence, the unit cell is doubled and the mirror 
plane is destroyed, resulting in the appearance of super- 



lattice peaks in the diffraction pattern. The new space 
group symmetry is i?3c, with hexagonal lattice param- 
eters aH ~ aoV2 and ch ~ 2aoV3. It should be em- 
phasized that these must be regarded as "average" long- 
range structures, since the presence of short-range or- 
der due to local displacements has been clearly demon- 
strated by the appearance of other types of superlattice 
peaks in electron diffraction studies2ii2iii not observed 
in x-ray or neutron diffraction patterns^i^ Significant 
deviations of the local atomic structure from the crys- 
tallographic long-range structure have also been found 
from pair-distribution function (PDF) analysis of time- 
of-flight neutron dataiiSiiSji* 

In the tetragonal region of the phase diagram, the 
space group is PAmm and the polar shifts lie along the 
[001] axis (flT w ct « ao,CT/aT > 1). Nevertheless, the 
time-of-flight neutron data show that this too should be 
viewed as an "average" long-range structure. In addition, 
Raman scattering studies have revealed the presence of 
local displacements of lower symmetry, which are also re- 
flected in a broadening of some of the x-ray diffraction 
peaksii^ii^ The nature of the local structure has been re- 
vealed in more detail from the PDF analysis described in 
Ref. 14, which shows that there are only gradual changes 
through the MPB, and suggests that the local environ- 
ment of each element remains relatively invariant of com- 
position. It is furthermore proposed that the population 
of local Pb displacements changes between the pseudocu- 
bic <100> and <110> directions as a function of the 
Ti/Zr ratio. This model is supported by recent theoreti- 
cal calculations in which the Pb distortions are identified 
as the determining factor for the average structure of the 
systempH 

In the original x-ray study by Noheda et al? the unit 
cell of the low-temperature monoclinic phase (now usu- 
ally designated Ma—) was found to be doubled with re- 
spect to the primitive cell, with the monoclinic a and h 
axes directed along the [110] and [110] axes of the lat- 
ter (om ~ ~ aoV^, ca/ ~ oq, space group Cm). 
Based upon the atomic positions determined from Ri- 
etveld analysis of the synchrotron x-ray data from PbZr 
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o.52Tio.4803, it was concluded^ that at 20 K the polar 
axis was tilted about 24° from the [001] axis towards the 
pseudocubic [111] axis. The structure can be regarded as 
a condensation of either the local displacements present 
in the tetragonal PAmm phase along one of the (110) 
directions, or alternatively those present in the rhombo- 
hedral i?3m phase along one of the (100) directions, as 
inferred by Corker et al.^ 

However, it is clear that there is a missing ingredient 
in this simple picture, for in a neutron powder diffrac- 
tion study of the same sample, Noheda et al. reported 
the presence of one very weak superlattice peak at 20 K 
corresponding to a doubling of the c axis of this mono- 
clinic cell, but did not identify the nature of this addi- 
tional distortionii& A similar cell-doubled phase was also 
observed for x — 0.48 by Ragini and coworkers in elec- 
tron diffraction patterns below 200 K, but not in their 
low-temperature x-ray patternSfl& Based on a subsequent 
Rietveld analysis of neutron powder data collected at 
10 K, the structure of this new phase was reported by 
Ranjan et al. to be monoclinic, with space group Pc}'^ 
The appearance of the weak superlattice reflections was 
attributed to antiphase tilting of the oxygen octahedra 
about the [001] direction, corresponding to an R-point 
instability in the cubic Brillouin zone. It was later re- 
ported that the correct space group for this proposed 
model was in fact CcjSS and a modified set of refined 
structural parameters was presented. 

The X = 0.48 composition has also been the subject 
of a recent low-temperature neutron powder study by 
Frantti and coUeaguesiSi They, too, note the presence 
of similar superlattice reflections, but reach very differ- 
ent conclusions; namely that these reflections are at- 
tributable to a minority rhombohedral phase with i?3c 
symmetry in coexistence with the monoclinic Cm phase, 
a model that was not considered by Ranjan et al}'^ or 
Hatch et alm^ In a footnote to their paper, Frantti et 
al. comment that the monoclinic Pc and Cc models pro- 
posed by the latter authors predict peaks that are not 
observed experimentally, and that the observed superlat- 
tice peaks can be accounted for by the R3c phase. How- 
ever, this conclusion was not supported by the results 
obtained by Noheda et alMi in an electron diffraction 
study of the same x = 0.48 sample used in the earlier x- 
ray study?^ which showed no evidence for a rhombohedral 
phase, but instead the monoclinic Cm phase in coexis- 
tence with nanoregions of a minority Cc phase ranging 
in size from 30-100 A. These conclusions have been ques- 
tioned by Frantti et al., who comment that their neutron 
data provide no evidence of a Cc phase, and argue that 
since electron diffraction probes only small volumes of 
the sample, it is generally not suitable for the determi- 
nation of average symmetry, and furthermore that the 
ion-milling technique used for sample thinning is a very 
violent one which can easily generate significant defects. 

In the light of these different interpretations, we have 
undertaken a detailed Rietveld analysis of the neutron 
data cited by Noheda et al& in an attempt to discrim- 



inate between the three models described above. Plau- 
sible results were obtained in all three cases, illustrat- 
ing how difficult it is to identify the correct structural 
model in complex systems of this type simply on the ba- 
sis of the standard goodness-of-fit criteria. Nevertheless, 
we conclude that, taken in conjunction with the electron 
diffraction data, the results point strongly towards the 
coexistence model of Cm and minority Cc phases. 



II. EXPERIMENTAL 

The sample consisted of about 4 g of sintered pellets 
roughly 1 cm in diameter and 1 mm thick from the same 
batch of material used in the previous x-ray studyi^ Long- 
range fluctuations in the composition of the x-ray sam- 
ple, Ax, were estimated to be less than ± 0.003 based 
upon an analysis of the peak widths. The pellets were 
loaded into a thin- walled vanadium can and mounted in a 
closed-cycle helium cryostat. Data were collected at the 
NIST reactor on the powder diffractometer BTl with a 
Cu monochromator set for a wavelength of 1.54 A, coUi- 
mation of 15' and 40' before and after the monochroma- 
tor, and 10' in front of each of the 32 ^He detectors. With 
this configuration, the best angular resolution attained is 
about 0.2° at 261 w 80° {Ad/d w 2 x IQ-^). 

Extended data sets were collected at 26 step intervals 
of 0.05° in the monoclinic region at 20 K, in the vicinity 
of the monoclinic-tetragonal transition at 325 K, and in 
the tetragonal region at 550 K. Analysis of the data was 
carried out with the FULLPROF program, ^■^ using the 
pseudo-Voigt peak-shape function with appropriate cor- 
rections for instrumental asymmetric broadening,^"* and 
linear interpolation between background points. Partic- 
ular attention was paid to the problem of anisotropic 
peak-broadening, which reflects the fact that closely- 
adjacent peaks may have markedly different widths aris- 
ing from local strains or compositional fluctuations, for 
example, as previously noted for PZT and related piezo- 
electric systems iSiSiiSiS^ In standard Rietveld analysis the 
peak widths are assumed to vary smoothly as a func- 
tion of scattering angle, and it is important to note 
that anisotropic peak-broadening due to microstructural 
effects can be mistakenly interpreted as a symmetry- 
lowering distortion of the unit cell of the average long- 
range structure. With the rapidly increasing use of high- 
resolution x-ray and neutron techniques, it is becoming 
clear that anisotropic peak-broadening is a common fea- 
ture of powder diffraction patterns, and should be al- 
lowed for as appropriate. One convenient way to do this 
is provided by the phenomenological model recently pro- 
posed by StephenS)S& in which the broadening is repre- 
sented by a series of coefficients J2hkl ^HKLh^ 
[H + K + L = 4), which take into acccount the Laue 
symmetry of the space group and are incorporated as re- 
finable parameters in the Rietveld program. For tetrag- 
onal Amm and monoclinic 2/m symmetry, there are re- 
spectively 4 and 9 such coefficients. 
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The data analysis is now described in detail for the 
tetragonal phase at 550 K, the monoclinic phase at 20 K, 
and the intermediate region at 325 K. 

A. 550 K 

All the peaks could be unambiguously indexed in terms 
of a tetragonal cell with a = 4.060, c = 4.100 A, except for 
two very weak peaks attributable to the vanadium sample 
holder, which were excluded from the analysis. Rietveld 
refinement was carried out with individual isotropic tem- 
perature factors assigned and the atoms placed in the fol- 
lowing positions of space group PAmm: Zr/Ti and 0(1) 
in 1(a) sites at 0,0, z; 0(2) in 2(c) sites at 0.5,0, z; and 
Pb statistically distributed among the 4(d) sites at x, x, 0. 
The Pb positions correspond to random displacements 
in the (110) directions away from the origin, as noted in 
the previous synchrotron x-ray study.^ The refinement 
converged rapidly and smoothly to a goodness-of-fit 
value of 1.30. However, as in the x-ray study, an exam- 
ination of the observed and calculated peak profiles re- 
vealed a number of systematic discrepancies indicative of 
anisotropic peak broadening, and additional refinements 
were carried out in which various combinations of the 
four possible coefficients were allowed to vary. A definite 
improvement was obtained when 5004 was varied alone 
{x^ = 1.20), but the results obtained with additional co- 
efficients were judged to be of dubious significance, and 
these coefficients were accordingly set to zero. This re- 
sult most likely reflects the sensitivity of the c lattice 
parameter to the presence of long-range compositional 
fluctuations. The flnal reflnement results based on this 
model are listed in Table I, and the proflle and difference 
plots are shown in Fig. 2. 



B. 20 K 

A series of refinements was carried out for each of 
the three models described in the Introduction, namely: 
single-phase Cc, two-phase Cm/Cc, and two-phase 
Cm/ R3c. However, instead of Cc, the non-conventional 
space group setting Jc was chosen, which has the distinct 
advantage of having a near-orthogonal unit cell closely re- 
lated to the Cm cell, in which the mirror plane is replaced 
by a c-glide plane and the c axis is doubled. The unit- 
cell axes are related via the transformation a/ = -cc, b/ 
= be, c/ = ac + Cc, where the subscripts refer to the 
unit cells of the Ic and Cc space groups respectively. In 
this setting, it is much easier to visualize the small dis- 
placements from the ideal Cm atomic positions. In the 
Ic cell, the Pb atom was chosen to lie at the origin, with 
Zr/Ti and three inequivalent O atoms in fourfold general 
positions at x, y, z and x, —y, 1/2 + z, and at the related 
body-center sites. The Zr/Ti and 0(1) atoms are in po- 
sitions similar to those in the Cm structure at x, 0, z, the 
main difference being that they are no longer required 



to lie on a mirror plane at ?/ = 0. The 0(2) and 0(3) 
atoms are in two sets of positions derived from the x,y,z 
sites and the symmetry-equivalent mirror plane sites at 
x,—y,z in the Cm structure. 



1. Single-phase Ic model 

In the first series of refinements the atoms were in- 
tially assigned the positions found in our previous x-ray 
studji^ with the exception of the 0(2) and 0(3) atoms, 
which were displaced from the ideal Cm positions by 
small shifts corresponding to antiphase tilting of the oxy- 
gen octahedra about the [001] axis, as assumed by Ran- 
jan et alMi The corresponding positions chosen for 0(2) 
and 0(3) in the Ic structure were x — S,y — S, z/2 and 
1/2 + X + 6,1/2 + y + 6,z/2, where 6 is the shift in the 
X and y directions due to tilting. It is important to 
note that with such a constrained-tilt, rigid-octahedron 
model, the x and y values assumed for the 0(2) and 0(3) 
positions (in this case the Cm values found in the pre- 
vious x-ray study) do not change in the course of the 
refinement. The resulting fit was reasonably good (x^ 
= 1.97), but inspection of the individual peak profiles 
once again revealed some significant discrepancies due to 
anisotropic peak broadening. As before, a distinct im- 
provement was obtained when the anisotropy coefficient 
5oo4 was refined (x^ = 1-78), but further refinements with 
various combinations of the other eight anisotropy coef- 
ficients gave only minimally improved fits, and the one- 
parameter anisotropy model was accordingly adopted for 
subsequent refinements. At this point, the constraints on 
the Zr/Ti and 0(1) y parameters were relaxed, but the 
shifts from the ideal positions at y = and the improve- 
ment in the overall fit were judged to be insignificant. A 
similar result was obtained when the Zr/Ti compositional 
parameter x was varied. 

Further refinements were performed, first with con- 
strained models corresponding to rigid-octahedron tilt- 
ing about the [111] and [110] axes respectively, and fi- 
nally, with all the constraints on the 0(2) and 0(3) po- 
sitions removed. In the latter case, the refinement pro- 
ceeded smoothly and converged rapidly to a set of po- 
sitions which were much closer to those of the [OOlj-tilt 
model than the other tilt models independent of which 
tilt model was used to provide the initial values of the 
positions. However, in none of these cases did the overall 
fit appear to be significantly improved, and we therefore 
conclude that the [001]-tilt model is a reasonable choice, 
although it is clearly not possible to rule out the other 
models on the basis of the present data. The final re- 
finement was therefore carried out for the constrained 
[001]-tilt model, but with the x and y values for 0(2) 
and 0(3) derived from the results for the unconstrained 
model. This refinement yielded a value of 1.52, with 
refined parameters as listed in Table II (column 1). Also 
shown are the values reported by Hatch et al^ trans- 
formed from Cc to Ic symmetry (column 3). From a 
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comparison of the two sets of atomic positions, it ap- 
pears that the constraints apphed by Hatch et al. do 
not in fact correspond to an [001]-tilt model, but instead 
to a simpler model in which only the y parameters of 
the 0(2) and 0(3) atoms arc displaced from their ideal 
Cm positions. An additional refinement based on such 
a y-shift model yielded results which are seen to be in 
excellent agreement with those of Hatch et al. (column 2 
of Table H) , although the fit is somewhat inferior to that 
given by the [001]-tilt model (x^ = 1.62). It is also worth 
noting that the values of axial ratio, Co/ag, and the pseu- 
docubic cell volume, Vq, obtained by the latter authors 
suggest a slightly higher Zr content (« 0.5%) relative to 
the present sampleiS^ 

From the atomic positions listed in the first column of 
Table II the octahedral tilt angle is calculated as about 
3°. The polar displacements of the Zr/Ti and Pb atoms 
with respect to the respective polyhedra centers are -0.08 
and -0.22 A along monochnic [100], and 0.18 and 0.44 A 
along [001], corresponding to a rotation of the polar axis 
towards pseudocubic [111] of roughly 25°. However, these 
values are representative only of the average long-range 
structure, since they do not allow for the local distortions 
revealed in the PDF analysis cited earlierii^ 

2. Two-phase Cm/R3c model 

The next set of refinements was performed for the 
two-phase Cm/R3c model favored by Frantti et alm^ 
Significantly better peak profiles were obtained with an 
anisotropic-broadening model for the Cm phase in which 
S'oo4 was allowed to vary, together with an isotropic 
particle-size broadening coefficient for the i?3c phase. 
The refinement converged rapidly to a value of 1.27 
with the final parameters as listed in Table HI. The lat- 
ter are in close agreement with those reported by Frantti 
et al., including the respective weight fractions of the 
two phases. Compared to the single-phase Ic model, the 
overall fit is considerably better {x^ = 1.27 vs. 1.52), 
but because several additional variable parameters are 
involved, it is difficult to judge the true significance of 
this result. Inspection of the results listed in Table III 
reveals that in both samples, Vq for the rhombohedral 
phase is larger by about 0.3 A'^, which would imply a 
significantly higher Zr content of some 3-4%iSiSi For the 
present sample, at least, such a conclusion would be in- 
consistent with the previously estimated long-range com- 
positional fiuctuationsi^ 

3. Two-phase C'm/Ic model 

The final set of refinements was carried out for the 
two-phase Cm/Ic model deduced by Noheda et al. from 
the results of an electron diffraction study.^^ Since we 
did not anticipate that a meaningful result would be 
obtained for an unconstrained refinement of two such 



closely-related structures, a highly constrained model 
was used; namely, the atomic positions in the Cm and 
Ic phases were constrained to be equivalent except for 
one additional parameter S for the latter representing 
the displacement along the x and y axes for the ideal- 
ized [001]-tilt model described above. The peak-shape 
model, including an 6*004 anisotropy coefficient, was also 
constrained to be equivalent for both phases, except for 
an isotropic particle-size broadening coefficient which was 
included for the Ic phase. The refined values for the two 
phases are listed side-by-side in Table IV for easy com- 
parison, and the profile fit and difference plot are shown 
in Fig. 3. The relative proportions of the Cm and Ic 
phases are approximately 4:1 and thus consistent with 
the electron diffraction results, but the estimated parti- 
cle size derived from the broadening coefficient is much 
larger, about 1000 A compared to 100 A. As pointed out 
by Frantti et al.,"^^ this discrepancy could arise because 
of the ion-milling techniques used to thin the electron 
diffraction sample, which can generate significant num- 
bers of defects. 

Detailed comparison of the results in Tables II, III and 
IV reveals that a better fit is obtained with the two-phase 
Cm/Ic model (x^ = 1.16, = 0.064) than with the 
Cm/R3c model (x^ = 1.27, i?^p = 0.067) or the single- 
phase Ic model (x^ = 1.52, Ru,p = 0.073), but it would 
nevertheless be premature to conclude that the former 
must therefore be correct, since there are no generally- 
accepted statistical tests to judge the true significance of 
the results. However, although the extended profile and 
difference plots for the two latter models are hardly dis- 
tinguishable by eye from those shown in Fig. 3, there are 
significant differences in some of the individual peak pro- 
files which provide additional insight, as shown in Figs. 
4 and 5 respectively. Fig. 4 shows the region around the 
strongest superlattice peak at 29 ft! 36.8°, from which it 
is evident that a much better fit is obtained with the Ic 
and Cm/Ic models than with the Cm/RZc model. On 
the other hand, the fit shown in Fig. 5 in the pseudocubic 
(200) region reveals serious deficiencies for the Ic model 
compared to the Cm/Ric and particularly the Cm/Ic 
model, which accounts much better for the asymmetry 
of the profiles in the central region. 

We note also that further analysis of the synchrotron 
x-ray data reported in Ref. 3 shows that the Cm/Ic co- 
existence model is superior to the two-phase Cm/ Pmim 
model previously used, with x^ values of 7.3 and 12.8 re- 
spectively. The results are in reasonable agreement with 
those of the neutron study; in particular, the ratio of the 
two phases is found to be about 4:1, very similar to the 
value listed in Table IV. The failure to detect any super- 
lattice peaks analogous to the one in the neutron pattern 
can be explained by the relatively much weaker x-ray 
scattering power of oxygen compared to Pb and Zr/Ti, 
resulting in calculated intensities that are insignificant 
compared to the background signal. 

From the atomic positions listed in Table IV, the oc- 
tahedral tilt angle about the [001] axis is calculated as 
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about 7°. The polar shifts of the Zr/Ti and Pb atoms 
with respect to the polyhedra centers are essentially the 
same as those obtained for the single-phase Ic refinement. 
It is also seen that the values of the cell volume Vq and 
the axial ratio cg/ao for the Ic phase are respectively 
slightly larger and smaller than those for Cm, and thus 
suggestive of a slightly higher Zr content 0.5%) for 
the former iSiSi 



C. 325 K 

Refinement was first carried out based on a model sim- 
ilar to that used for the 550 K data; namely, a single- 
phase tetragonal structure with PAmm symmetry, Pb 
atoms statistically distributed among the 4(d) sites at 
x,x,0, and a single anisotropy-broadening coefficient, 
S'oo4- However, the overall fit was only mediocre ( 
— 2.94, Rujp — 0.072), and a detailed inspection of the 
individual peak profiles revealed asymmetries consistent 
with the presence of a monoclinic component. Such a co- 
existence model of monoclinic and tetragonal phases for 
X = 0.48 at room temperature was proposed in an earlier 
neutron study by Frantti et al.^ and in a more recent 
x-ray study by Ragini et ali^ Further refinements based 
on this model gave a markedly improved fit ( = 1.76, 
Rwp = 0.056), but some residual diffuse scattering was 
clearly present between some of the peaks. This scatter- 
ing is probably associated with locally disordered regions 
in the vicinity of domain walls and can be modeled in 
a simple, albeit rather artificial, way by the addition of 
a cubic phase with Pm3m symmetry, as assumed in our 
previous x-ray studyi^ Such a three-phase model yielded 
a reasonably satisfactory fit (x^ = 1.47, R^p = 0.051), 
with weight fractions of tetragonal, monoclinic and cubic 
phases in the ratio 0.61:0.33:0.06. 

The refined parameters are listed in Table V, and the 
profile fit and difference plot are shown in Fig. 6. Also 
listed are the parameters reported by Ragini et ali^ and 
Frantti et alm^ In the latter case, it is rather surprising in 
the light of the results reported in Ref. 6 that the lattice 
strain co/ap was found to be significantly larger for the 
monoclinic phase than for the tetragonal one, since one 
would not expect rotation of the polarization direction 
away from [001] in the monoclinic phase to increase this 
strain. Other than this, the three sets of parameters are 
in reasonable agreement except that the fraction of Cm 
phase in the room temperature studies is considerably 
larger than at 325 K, as would be expected. 

Further analysis of the earlier x-ray dataSi revealed that 
this three-phase model gives a noticeably better profile 
fit than that obtained with the two-phase PAmm/ Pmim 
model previously used, with values of 7.5 and 9.9 
respectively. The weight fractions of the three phases 
were in the ratio 0.55:0.40:0.05, comparable to the neu- 
tron values listed in Table V. We emphasize, however, 
that these results should be regarded as representative 
only of an average long-range structure, since the true 



nature of the material in the transition region is surely 
far more complex than implied by a simple three-phase 
model. It is more likely in this temperature interval that 
small fluctuations in composition lead to coexistence of 
the tetragonal phase with locally-ordered monoclinic re- 
gions of widely varying sizes and possibly some disor- 
dered regions. 



III. DISCUSSION 

In summary, the results obtained in the present 
neutron investigation are consistent with the coex- 
istence of majority Cm and minority Ic phases in 
PbZro.52Tio.48O3, in agreement with the results of a re- 
cent electron diffraction study of the same sample. The 
Cm (Ma) phase, which plays a key role in the piezo- 
electric and ferroelectric behavior of PZT and related 
systems is the majority phase at low temperature. The 
structure of the minority Ic phase is readily visualized 
as the superposition of an antiphase octahedral-tilt sys- 
tem on the parent Cm structure. Furthermore, the close 
agreement between the refinement results for the alter- 
native Ic and Cm/Ric models and those in Refs. 20 
and 21 respectively suggests that the two-phase Cmjlc 
model is worth consideration in those cases as well. It is 
possible that the coexistence of Cm and Ic phases in the 
present x = 0.48 sample reflects the existence of a nar- 
row thermodynamically-stable region with Ic symmetry 
at low temperature somewhere between 0.45< x < 0.48. 
In this case, the coexistence of Ic and Cm phases could 
plausibly be attributed to the presence of long-range 
compositional fluctuations, as suggested by the values 
of the lattice parameters for the two-phase refinement in 
Table V. In this context, it is interesting to note that re- 
cent neutron data obtained by Frantti et al. for a sample 
with X — 0.46 show clear evidence of a superlattice peak 
at 4 K.^^ However, the authors interpret this as evidence 
for the coexistence of Cm and i?3c phases, and did not 
consider the possibility of Cm and Ic phase coexistence. 
Alternatively, the Ic phase could be a metastable one re- 
sulting from the presence of local strains at domain-wall 
boundaries, for example. Indeed, it is noteworthy that 
first-principles calculations by Fornari and Singh have 
shown that local stress fields may lead to the coexistence 
of both ferroelectric and rotational instabilities near the 
MPB.-^^ In any case, it is clear that a very careful high- 
resolution x-ray, neutron and electron diffraction study of 
extremely well-characterized samples would be required 
in order to throw further light on these issues. 

On a final note of caution, the present study also 
demonstrates that the interpretation of the results of Ri- 
etveld analysis may be very tricky for complex systems 
such as this one in which allowance must be made for 
the possible coexistence of closely-related phases and the 
presence of anisotropic peak broadening. The choice of 
any particular model should take into account not only 
the quality of the refinement as judged by the agreement 
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factors and goodness-of-fit, but also the diffraction pro- 
files of alternative models in selected key regions of the 
pattern, and, if feasible, data from complimentary struc- 
tural techniques such as electron diffraction. 
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TABLE L Refined structural parameters for tetragonal 
PbZro.52Tio.48O3 at 550 K, space group PAmm, lattice pa- 
rameters a = 4.0596(1), c = 4.0999(1) A. The refinement was 
based on a model with the Pb atoms statistically distributed 
among 4(d) sites at x,x,Q, corresponding to local displace- 
ments along (110) directions. Figures in parentheses denote 
standard errors referred to the least significant di git. Rwp, 
Rb and are agreement factors as defined in Ref. l2a 



X y 
Tb 0.033(1) 0.033(1) 
Zr/Ti 0.5 0.5 
0(1) 0.5 0.5 
0(2) 0.5 0.0 



0.0 0.028(1) 
0.450(2) 0.005(1) 
-0.061(1) 0.027(1) 
0.427(1) 0.027(1) 



Rwp 

Rb 

X2 



0.048 
0.034 
1.20 



TABLE IL Refined structural parameters for monoclinic 
PbZro.52Tio.48O3 at 20 K, single-phase model with space 
group /c, for the [001]-tilt and y-shift models described in 
the text. The Pb atom was fixed at the origin, and the 0(2) 
and 0(3) temperature factors were constrained to be equal. 
Also listed are the parameters recently reported by Hatch 
et alw~, but with the values transformed from Cc to Ic sym- 
metry. Vo and co/0,0 represent respectively the volume and 
axial ratio of the primitive pseudocubic cell, with cq — c/2 
and ao ^ {a + b)/2^2. 



Present study Hatch et al. 







[001] -tih 


y— shift 


y— shift 


a(A) 




5.7131(1) 


5.7131(1) 


5.7312(7) 


6(A) 




5.7000(1) 5.7001(1) 


5.7093(6) 


c(A) 




8.2679(2) 


8.2683(3) 


8.2363(7) 






90.475(2) 90.473(2) 


90.50(1) 




) 


67.31 


67.31 


67.37 


Col Ho 




1.0245 


1.0246 


1.0181 


Pb: 


c/(A2) 


0.013(1) 


0.012(1) 


0.013(1) 


Zr/Ti 


: X 


0.524(2) 


0.524(2) 


0.519(5) 






0.219(1) 


0.218(1) 


0.216(2) 






0.002(1) 


0.003(2) 


0.006(4) 


0(1): 


X 


0.542(1) 


0.543(1) 


0.548(3) 






-0.046(1) 


-0.046(1) 


-0.044(1) 






0.011(1) 


0.011(1) 


0.011(3) 


0(2): 


X 


0.275(1) 


0.287(1) 


0.289(2) 




V 


0.243(1) 


0.233(1) 


0.233(1) 






0.193(1) 


0.194(1) 


0.196(1) 






0.010(1) 


0.011(1) 


0.009(1) 


0(3): 


X 


0.801(1) 


0.787(1) 


0.789(2) 




y 


0.768(1) 


0.767(1) 


0.767(1) 






0.193(1) 


0.194(1) 


0.196(1) 






0.010(1) 


0.011(1) 


0.009(1) 






0.073 


0.076 


0.086 


Rb 




0.041 


0.047 


0.040 


x' 




1.52 


1.62 


1.21 



TABLE III: Refined structural parameters for 
PbZro.52Tio.48O3 at 20 K, two-phase model with space 
groups Cm and Ric. The Pb and Zr/Ti atoms were fixed 
at the origin for the Cm and J?3c refinements respectively, 
and the temperature factors for the separate atoms were 
constrained to be the same in both structures. Also listed 
are the 10 K parameters recently reported by Frantti et aim— 







Present study 


Frantii et al. 






Cm 


R3c 


Cm 


R3c 


a{ A) 




5.7120(1) 


5.7415(6) 5.7097(7) 


5.744(2) 


b( A) 




5.6988(1) 




5.6984(7) 




c( A) 




4.1353(1) 


14.208(3) 4.1367(3) 14.212(8) 


pn 




90.479(2) 




90.449(8) 






67.32 


67.60 


67.29 


67.68 


Co/ao 




1.0257 


1.0 


1.0256 


1.0 


Pb: 






0.283(4) 




0.282(5) 






0.012(1) 


0.012(1) 


0.009(1) 


0.004(6) 


Zr/Ti: 


X 


0.531(2) 




0.539(3) 








0.441(2) 




0.441(3) 






u{K') 


0.004(1) 


0.004(1) 


0.001(2) 


0.001(2) 


0(1): 


X 


0.543(1) 




0.540(1) 








-0.090(1) 




-0.092(2) 






[/(A^) 


0.008(1) 




0.011(2) 




0(2): 


X 


0.288(1) 


0.137(3) 


0.283(1) 


0.148(3) 




y 


0.254(1) 


0.347(3) 


0.253(1) 


0.354(3) 






0.389(10 


0.081(4) 


0.388(1) 


0.081(6) 






0.015(1) 


0.015(1) 


0.013(1) 


0.007(5) 


/(wt fraction) 


0.89(1) 


0.11(1) 


0.87 


0.13 






0.067 


0.058 


Rb 




0.039 


0.067 


0.042 






1.27 


2.28 



TABLE IV: Refined structural parameters for monoclinic 
PbZro.B2Tio.48O3 at 20 K, two-phase model with space groups 
Cm and 7c. The atomic positions for 7c symmetry were based 
upon the [001]-tilt model described in text and constrained 
to be equivalent to those for Cm except for one additional 
parameter 5 corresponding to 0(2) and 0(3) displacements 
in the x and y directions due to tilting. For comparison 
with the 7c structure, the 0(2) and 0(3) positions for the 
Cm. structure arc shown separately, although in fact they are 
synunetry-equivalent. The temperature factors for the sepa- 
rate atoms were constrained to be the same in both structures. 







Cm 


7c 


«(A) 




5.7097(1) 


5.7401(7) 


&(A) 




5.6988(1) 


5.7188(8) 


c(A) 




4.1373(1) 


8.2098(11) 






90.473(2) 90.550(10) 


K(A3) 


67.31 


67.37 


Cojao 




1.0257 


1.0127 


Pb: 




0.012(1) 


0.012(1) 


Zr/Ti: 


X 


0.530(2) 


0.530(2)- 






0.437(1) 


0.218(1) 




f/(A2) 


0.003(1) 


0.003(1) 


0(1): 


X 


0.541(1) 


0.541(1) 






-0.089(1) 


-0.045(1) 






0.011(1) 


0.011(1) 


0(2): 


X 


0.286(1) 


0.257(1) 




y 


0.254(1) 


0.225(1) 






0.390(1) 


0.195(1) 




c/(A2) 


0.015(1) 


0.015(1) 


0(3): 


X 


0.786(1) 


0.814(1) 




V 


0.754(1) 


0.783(1) 






0.390(1) 


0.195(1) 






0.015(1) 


0.015(1) 


/(wt fraction) 


0.78(2) 


0.22(2) 






0.064 


Rb 




0.039 


0.053 






1.16 
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TABLE V: Refined structural parameters for the Cm and 
PAmm phases in PbZro.52Tio.48O3 at 325 K with the three- 
phase model described in text. The temperature factors for 
the separate atoms were constrained to be the same in both 
structures. Also listed are the room-temperature parameters 
recently reported in an x-ray study by Ragini et alw~ and in 
a neutron study by Frantti et al^ The weight fraction of the 
cubic Pm3m phase was determined as 0.06. 





Present study 


Ragini 


et al. 


Frantti et al. 




Cm 


PAmm 


Cm 


PAmm 


Cm 


PAmm 


a(A) 


5.7268(3) 


4.0393(1) 


5.7520(1) 4.0429(2) 


5.7129(3) 


4.0550(4) 


6(A) 


5.7187(3) 


- 


5.7431(2) 




5.7073(3) 


- 


c(A) 


4.1230(2) 


4.1388(1) 


4.0912(4) 4.1318(3) 


4.1436(1) 


4.1097(6) 


pn 


90.393(5) 




90.48(1) 




90.199(3) 






67.51 


67.53 


67.57 


67.53 


67.55 


67.58 




1.0189 


1.0246 


1.0067 


1.0219 


1.0262 


1.0135 


Pb: z 




0.035(2)" 






- 


- 




0.017(1) 


0.017(1) 


0.107'' 


0.030'= 


0.02T^ 


0.019(1) 


Zr/Ti: x 


0.530(4) 




0.578(3) 




U.oO / (2 j/0. 494(4 ) 






0.432(5) 


0.442(2) 


0.473(3) 


0.447(2) 


0.426(l)/0.404(4)" 


0.431(4) 




0.003(1) 


0.0043(1) 


0.015(1) 


0.005(2) 


0.004(1) 


0.019(1) 


0(1): X 


0.540(2) 




0.50(1) 




0.522(1) 






-0.080(3) 


-0.085(2) 


-0.10(1) 


-0.109(6) 


-0.090(1) 


-0.080(2) 




0.016(1) 


0.016(1) 


0.00(1) 


0.029(1) 


0.013(1) 


0.019(1) 


0{2) : X 


0.287(2) 




0.36(1) 




0.270(1) 




y 


0.255(1) 




0.219(8) 




0.252(1) 






0.400(2) 


0.0395(1) 


0.404(8) 


0.389(3) 


0.391(1) 


0.400(1) 




0.021(2) 


0.021(2) 


0.04(1) 


0.029(1) 


0.013(1) 


0.019(1) 


f (wt fraction) 


0.33(1) 


0.61(2) 


0.58 


0.42 


0.69 


0.31 


Rwp 


0.051 


0.128 


0.021 




Rb 


0.043 


0.030 


0.041 


0.062 








1.47 


3.39 


2.69 





(a) Pb atoms statistically distributed among 4(d) sites at x, 0. 

(b) Equivalent isotropic U 

{Uii = 0.221, U22 = 0.027, (733 = 0.074, U13 = 0.030 A^). 

(c) Equivalent isotropic U (C/11 = U22 = 0.031, U33 = 0.027 A^). 

(d) Equivalent isotropic U 

(C/11 = 0.027, U22 = 0.026, U33 = 0.011, U13 = 0.013 A^). 

(e) Zr and Ti parameters refined independently. 
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FIG. 1: PZT phase diagram as originally proposed by JaflFe 
et al. in Ref. [l| (open circles) with the modifications reported 
by Noheda et al. in Ref. ^ (full circles). The various phases 
described in the text are denoted by C (cubic Pm3m), Rht 
(rhombohedral Rim), Rlt (rhombohedral _R3c), T (tetrago- 
nal PAmm), and Ma (monoclinic Cm). 
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FIG. 2: (Color) Observed and calculated diffraction profifes 
from the Rietveld refinement of PbZro.52Tio.48O3 at 550 K, 
with space group PAmm. The difference plot is shown be- 
low, with short vertical markers denoting the calculated peak 
positions. 
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FIG. 3: (Color) Observed and calculated diffraction profiles 
from the two-phase Rietveld refinement of PbZro.52Tio.48O3 
at 20 K, with space groups Cm and Ic. The difference plot 
is shown below, with upper and lower sets of vertical mark- 
ers denoting the calculated peak positions for Cm and Ic 
respectively. The position of the weak superlattice peak at 
28 ^ 36.8° (pseudocubic 3/2 1/2 1/2) is indicated with an 
asterisk. 



20K 




36 37 38 39 36 37 38 39 36 37 38 39 40 



2Q^deg) 20(de.) 2e(deg) 



FIG. 4: (Color) Observed and calculated difltraction profiles 
and difference plots in the region around the strongest super- 
lattice peak from PbZro.52Tio.48O3 at 20 K for single-phase 
Ic (left panel), two-phase Cm + R3c (center panel), and two- 
phase Cm + Ic (right panel) . 
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FIG. 5: (Color) Observed and calculated difltraction profiles 
and diflterence plots in the region around the pseudocubic 
(200) reflection from PbZro.52Tio.48O3 at 20 K for single- 
phase Ic (left panel), two-phase Cm + RSc (center panel), 
and two-phase Cm + Ic (right panel) . 
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FIG. 6: (Color) Observed and calculated diffraction profiles 
from the three-phase Rietveld refinement of PbZro.52Tio.48O3 
at 325 K, with space groups PAmm, Cm and Pm3m. The 
difference plot is shown below, with upper, middle and lower 
sets of vertical markers denoting the calculated peak positions 
for PArara, Cm and Pm3m respectively. 



